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Although nonclassical quantum states are important both conceptually and as a resource for 
quantum technology, it is often difficult to test whether a given quantum system displays nonclassi- 
cality. A simple method to certify nonclassicality is introduced, based on easily accessible collective 
atomic quadrature measurements, without the need of full state tomography. The statistics is an- 
alyzed beyond the ground-state noise level, by direct sampling of a regularized atomic quadrature 
quasiprobability. Nonclassicality of a squeezed ensemble of 2 • 10 Cesium atoms is demonstrated, 
with a significance of up to 23 standard deviations. 
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Introduction. The difference between the classical 
and quantum mechanical description of physical systems 
is often discussed in terms of nonclassicality. This topic is 
of general interest, since it gives insight into the reasons 
why a classical picture of nature has to fail, and provides 
methods to explore what a quantum mechanical system 
can do which a classical system cannot. In particular, 
such issues are of great importance for numerous fields 
in quantum optics, quantum computation, quantum in- 
formation, and related research areas. 

Nonclassicality of the harmonic oscillator has been well 
understood during the last years. Its definition is based 
on the coherent states \a), parameterized by a complex 
amplitude a. The evolution of these pure quantum states 
most closely corresponds to the classical evolution of po- 
sition and momentum. They are minimum uncertainty 
states with a stable Gaussian wave packet. In this sense 
they represent the closest analogue to the classical behav- 
ior. Furthermore, they allow to decompose any quantum 
state in the form 

p = J d 2 aP(a)\a)(a\, (1) 

where P(a) is the Glauber-Sudarshan P function [TJ[2]. 
In this representation, an arbitrary quantum state for- 
mally resembles a statistical mixture of (classical) coher- 
ent states. However, for many states P(a) may show neg- 
ativities, which cannot be interpreted in terms of statis- 
tical probability densities. Such states are consequently 
referred to as nonclassical [3]. 

Atomic coherent states (ACS) have already been de- 
veloped in [4]. It has been shown that their interaction 
with a classical electric and magnetic field can be un- 
derstood classically, cf . [5 , and any density operator can 
be formally written as a statistical average over ACS. 
This remarkable property allows one to define nonclassi- 
cality in a way closely corresponding to that of the har- 
monic oscillator [6] [7] . Specific nonclassical properties of 
atomic systems have already been discussed in previous 



publications, such as spin squeezing [5HTU] and squeezing 
in Ramsey spectroscopy [TT]. The creation of entangle- 
ment in spin systems |121 1 13] and the relation between 
entanglement and squeezing [14] have also been studied. 
Recently, nonclassicality has been examined in terms of 
measurable probability distributions 15J. 

The present letter addresses a method to visualize 
nonclassicality of atomic systems, characterized by an- 
gular momentum states with a large orientation in the 
z-direction and small fluctuations along the two orthog- 
onal directions. To provide a simple and experimentally 
accessible condition for the detection of the nonclassi- 
cal features of the system, we introduce the concept of a 
regularized atomic quadrature quasiprobability (AQQP). 
Our method is applied to experimental data from ap- 
proximately 10 5 two-level atoms, which were collectively 
prepared in an atomic spin-squeezed state. It uncovers 
the nonclassical properties, even for weak squeezing. 

Atomic coherent states. Formally, the atomic ensem- 
ble under study can be characterized as an angular mo- 
mentum system. In quantum mechanics, the angular 
momentum is described by a vector valued operator 

J = (Jx, Jy, Jz), satisfying the commutation relations 
(setting h = 1): 

[Jx: Jy] — '^J z, [Jy, Jz] — iJx, [Jz, Jx] ^Jy (^) 

Furthermore, we introduce the angular momentum lad- 
der operators J± = J x ±iJ y , whose commutator is given 
by j_] = 2J Z . Throughout this work, we keep the 
quantum number j fixed, while the quantum number m 
takes values in {—j, — j + 1, . . . 

The foundation for the discussion of nonclassicality 
is typically given by the coherent states, being defined 
as pi] 

\B,<p) = exp{(^J_ - 6e-^J+)/2} , (3) 

where the fiducial state \j',j) is the common eigenstate 
of the operators J 2 and J z with maximum eigenvalue for 



2 



the latter operator. These coherent states can be seen as 
the analog to a classical angular momentum state. They 
have minimum uncertainty in J 2 , and their expectation 

value of the projection of J, 

n(9, <p) ■ J = sin(0) cos((/j) J x + s'm(9) sm(ip)J y + cos(9)J z 

(4) 

is exactly equal to j, without uncertainty, such that the 
state \9, (p) can be understood as pointing exactly in the 
direction of n(9,tp). Moreover, their time evolution in 
a classical magnetic field exactly matches the classical 
evolution [5] . Finally, there exists a representation of the 
density operator of an arbitrary quantum state in the 
form 

P = £ sin Od0 JJ P(9, <p) \9, v ) (9, <p\ . (5) 

Similar to Eq. ([I]), this is formally equivalent to a decom- 
position of a quantum state into a statistical mixture of 
(classical) coherent states. Therefore, all quantum states, 
for which the weight function P(9, if) is non-negative, 
are classical angular momentum states. Conversely, any 
quantum state for which such a non-negative function 
P(9,ip) does not exist, will be referred as being nonclas- 
sical 017]. 

The Holstein-Primakoff-approximation. Let us con- 
sider a system with a large angular momentum, i.e. j 3> 
1, which is mainly aligned in the z-direction. In this 
case, the operator J z can be approximated by its expec- 
tation value ( J z ) ss j, being a real number as in classical 
physics. Only the remaining angular momentum compo- 
nents J x , J y are treated fully quantum-mechanically [IT] . 
Then, the ladder operators satisfy the approximated 
commutator relation 

[J+,J-]«2j. (6) 

The re-normalized operators j± — -^jJ± obey the 
same algebra as the bosonic operators. We may introduce 
rescaled angular momentum components via j± = g (j x ± 
ij y ), which are defined such that the uncertainty of the 
ground state, (Aj x ) 2 ) \j;j) = (j; j\ {Aj y ) 2 ) \j; j) = 
1. The weight function P(9,(p) corresponds to the 
Glauber-Sudarshan P function of the quantum state, 
with the complex argument a = yj j /29e lip . The proper- 
ties of P(a) form the basis for defining nonclassicality of 
the harmonic oscillator. Consequently, we can reformu- 
late established methods for nonclassicality for angular 
momentum systems. 

Atomic nonclassicality quasiprobabilities. Although 
the P function uniquely defines the nonclassicality of a 
quantum system, its reconstruction is in general impos- 
sible for the harmonic oscillator, since it can be highly 
singular. For getting definite answers about nonclassical 
effects, one may introduce a filtered atomic P function, in 



close analogy to that of the harmonic oscillator [IS] . For 
atomic systems, one multiplies the characteristic function 
$(0 dlP{B,<p), 

$(0 = (e^-~«*J+)e l « |2/2 , (7) 

with a suitable filter function fi^(^) (parameterized by 
w), 

* n (0 = *(e)n«.(0- ( 8 ) 

The atomic nonclassicality quasiprobability is then given 
by 

PaUxJy) = \ I dV W ' Cp - J " W 3n(0, (9) 
^" Jc 

with £ = £ r + The filter function f2„;(£) has to sat- 
isfy certain properties: First, it has to guarantee that 
the nonclassicality quasiprobability of any state is always 
regular. Second, it only should show negativities if the 
quantum state is nonclassical. For the harmonic oscil- 
lator we have shown (TS] that for any nonclassical state 
one finds a nonclassicality quasiprobability demonstrat- 
ing negativities by simply varying a real parameter w 
that scales the filter width. 

Atomic quadrature quasiprobability. In the following 
we will introduce the AQQP to efficiently characterize 
nonclassicality of atomic spin systems with a minimum 
of measured data. Quadrature distributions are just 
probability densities, so that a regularization procedure 
is superfluous. The situation changes if one analyzes 
the quadrature statistics beyond the ground-state noise 
level [T5], to visualize usually hidden nonclassical effects, 
also beyond the negativities of the Wigner function |20j . 
Subsequent filtering yields the marginals of Pn(jx,jy)- 

Introducing atomic quadrature operators j v , 

U = U^+i-e-' lv , (10) 
we define the AQQP via 
1 f°° 

p n ( Jip ) = —j dA$(Ae^)rUA)e iA ^, (11) 

which is just a marginal of Pn(j x ,j y ) for a fixed angle 
ip. Clearly, if such a marginal distribution shows nega- 
tivities, then Paijxijy) has negativities as well, and the 
state is nonclassical. In general, the inverse conclusion is 
not true. 

The AQQP distribution can be sampled from a set of 
N measured angular momentum values {jk)k=i directly 
by using a similar approach to [5T]. Hence an estimator 
for the AQQP is expressed as an empirical mean of a 
suitable pattern function fn(J;j v ,w), 

1 N 

Pn(j v ) = X] fnUk;3>p,w). (12) 
fe=i 
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The required pattern function is given by 

1 r°° 

fa(j',3 v ,w) = ir e k2 / 2 e ik( i-^n w (k)dk. (13) 

^ J-oo 

Furthermore, the uncertainty of pn(jip) can be esti- 
mated as l/V-ZV times the empirical standard deviation 
of the sampled numbers {/nC/fc! jtp, w)}^ =1 . Therefore, 
the sampling method is a very simple way to estimate 
the quantity of interest together with its statistical error. 

Experimental procedure. The experimental setup is 
described in detail in [23]: A cloud of N a w 2 ■ 10 5 
cold Cesium atoms is prepared by laser cooling in a 
magneto-optical trap (MOT) and subsequently trans- 
ferred into a far off-resonant dipole trap formed in the 
40 /im waist of a single 5W, 1064 nm laser beam. Em- 
ploying a sequence of optical pumping pulses, microwave 
pulses, (optical) purification pulses, a magnetic bias field 
of approximately 1.5 Gauss, and a resonant microwave 
7r/2 pulse, each atom is prepared in the state \iJ)q) — 
(\F=3,m F =0) + \F=4,m F =0))/V5. 

In a simplified picture we describe each atom as a two- 
level system with pseudo-spin 1/2 and with basis states 
||, ±|) := ^(\F = 3)±e iu " t \F = 4)), with frio 3i = 
A-Ehyperfmc- The atomic population difference in the two 
hyper-fine manifolds F = 3 and F = 4 now can be re- 
garded as the a;— component of a collective j = N a /2- 
pscudospin vector. This corresponds to the sum of the 
individual pseudo-spins, \j;j) = <8)^i |^o)»j an ACS. 

By detecting the differential phase shifts <\> that the 
atomic cloud imprints on a bichromatic, off-resonant op- 
tical probe pulse [23] we measure the atomic popula- 
tion difference AN between the two hyper-fine levels. 
In this way we perform a non-destructive quantum non- 
demolition (QND) measurement of j x , which projects the 
atomic state into a squeezed state around C4*i- Here 
4>i denotes the outcome of the QND measurement and 
the parameter £ is a function of the QND measurement 
strength, which is directly proportional to the optical 
density of the ensemble and the probability of sponta- 
neous photon emission (22| . A subsequent verification 
measurement <p2, performed on the same ensemble, shows 
sub-projection noise quantum fluctuations of the differ- 
ential phase shift with respect to the prediction ^(j>\. 

As it is characteristic for QND measurements the 
signal-to-noise ratio (and therefore the measurement 
strength) is limited by the quantum noise of the meter 
system - in our case by the shot noise 6n of the probe 
light: 

4> x = — + n&N, 2 = — + nAN, (14) 
ni n 2 

where k is the light-atom coupling strength. Although 
the destructive second measurement uses more photons 
(ri2 = 1.5 ni) than the entangling QND measurement, it 
still contains a fair amount of light shot noise (see FigfTI) . 
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Figure 1. Calibration procedure: Blue points: var(0 ACS ) for 
different atom numbers; each data point is formed by 4841 in- 
dividual measurements. The scaling of the noise as a function 
of the number of atoms N a allows us to identify noise con- 
tributions. Blue area: light shot noise. Green area: atomic 
projection noise. Red areas: technical fluctuations. Black 
dashed line: predicted noise for a ACS with the same collec- 
tive j z as our squeezed state. Magenta diamonds: var(0i2) 
for experiments with different N a . 

Data analysis. To be able to convert the N measure- 
ment outcomes 4>i2k '■— 4>2k ~ C 4>ik (k = X, . . . , N) into 
quadrature samples, we calibrate our experiment as fol- 
lows: 

• We prepare a ACS using a varying number of 
trapped atoms N a and, for calibration purposes 
only, subtract measurement outcomes <^ cs from 
successive independent MOT loading cycles (~ 5 s 
apart) to compensate for slow drifts of experimental 
parameters. By performing a scaling analysis [2"2"] 
on the differential values </> ACS = (^k° s —(pk+i )A/2 
we can identify the amount of light shot noise 
\&x(Sn/n) and atomic quantum projection noise 
var(/cAiV) and discriminate those against technical 
noise sources. 

• Additionally, due to spontaneous emission pro- 
cesses caused by the (/>i-measurement, performed 
with ni — 4.1 • 10 7 photons, the atomic state is 
not pure anymore and compared to the ACS the 
Ramsey-fringe contrast (and thus the macroscopic 
spin j of the ensemble) is reduced by a factor 
i] = exp(-me) = 65.8%, where e = 1.02 ■ 10" 8 
is the decoherence per photon, determined experi- 
mentally as described in [2"3"] . 

With these calibrations we can calculate the phase fluc- 
tuations that we would measure in the absence of any 
technical noise for a ACS comprised of ?7iV a atoms, which 
yields the same Ramsey fringe amplitude as our squeezed 
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state: 



var ACS ( J)A r a ) = var(fc 2 /n 2 ) + K 2 r)N a . 



(15) 



For a varying number of atoms < N a < 2.9 • 10 5 
we perform in total N — 84730 experiments and obtain 
phase shift measurements 012^, now without taking the 
difference between subsequent MOT loading cycles. 

We convert the 4>i2k m to quadrature values by normal- 
izing to va,T AC s( v N a y- 



3k = 012fc/V var ACS(r)JV a )- 



(16) 



The jk do not only contain the atomic signal jk but also 
added photon shot noise, which is Gaussian noise, uncor- 
rected with the atomic population difference AN. Due 
to the normalization procedure the effect of such added 
Gaussian noise is equivalent to a detection with a reduced 
quantum efficiency |24j (mixing with an ACS). The effec- 
tive quantum efficiency of our quadrature measurement 
depends on the 'atomic signal to shot noise'-ratio and is 
therefore a function of N a : 



Efficiency 



var ACS(r)Af a ) — v arACS(0) 

var ACS ( l;A , a ) 



(17) 



In our experiment, it reaches values up to 83% for mea- 
surements taken with the highest number of atoms N a = 
2.9 • 10 5 . This effect only attenuates negativities of the 
nonclassicality quasiprobability Pq and any nonclassical- 
ity found in jk is therefore sufficient to prove a nonclas- 
sical atomic state in the pseudospin operators j x , j y . 

Verification of nonclassicality. In the following, we 
only consider experiments performed with an effective 
detection efficiency above 77%, which show a quadra- 
ture variance 1.67 dB below the projection noise. From 
the {jfc}^] 4841 data points, we can estimate the AQQP 



according to Eq. (12), by using a one-dimensional auto- 
correlation filter 

Q w (k) =Af~ 1 I uj(k')uj(k' + k/w)dk, (18) 



with cj(fc) = e~ k and N = J[uj(k)] 2 dk. This filter pre- 
serves all information about the distribution, but only 
displays negativities for nonclassical states. First, we 
look at the significance of the negativity as a function 
of the filter width, 



E(u,) 



o-(pnOV)) 



(19) 



where o~(pn(j v )) is the standard deviation of the esti- 
mate pn(jtp)- This quantity is nothing but the value of 
the negativity in terms of standard deviations of the mea- 
surement. For observing nonclassical effects, this quan- 
tity has to be negative. From our experimental data, we 
always observe such negative significance. Moreover, the 
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Figure 2. AQQP distribution of a squeezed angular momen- 
tum state for a filter width w — 1.1. The error bars are 
smaller than the linewidths. The dashed green curve shows 
the standard quadrature distribution (i.e. the marginal of the 
Wigner function), which is non-negative. Inset: Significance 
E of the negativity of the AQQP versus the width parameter. 



inset in Fig. [2] shows the dependence of its absolute value 
on the filter width w. It can be seen that there is a filter 
width, for which the significance is optimized. In case of 
squeezing, this may even be arbitrarily close to w = 0. 
Here, we can obtain a significance of up to 23 standard 
deviations. 

The main part of Fig. [2] shows the AQQP of the cor- 
responding state. Here, the filter width is chosen to be 
w = 1.1, which gives a smaller significance of 10.0 stan- 
dard deviations, but also leads to more pronounced ef- 
fects. We observe oscillations with significant negativi- 



ties. Since the autocorrelation function (Eq. (18)) fulfills 
the requirements of a filter function as listed in the the- 
ory part by construction, the negativities of the AQQP 
clearly reflect the nonclassicality of the quantum state. 

Conclusions. We introduced the atomic quadrature 
quasiprobability, pn(j v ), for atomic systems described 
by angular momentum states in the Holstein-Primakoff 
approximation. Its negativities certify negativities of the 
Glauber-Sudarshan P-function of the atomic system, and 
hence the nonclassicality of the state beyond both the 
ground-state noise level and negativities of the Wigner 
function. We apply our method to experimental data ob- 
served in a squeezed atomic ensemble and demonstrate 
nonclassicality with a statistical significance of up to 23 
standard deviations. We consider our method a valuable 
tool for certifying nonclassicality especially for states pre- 
pared in composed many-identical-particle systems such 
as cold atomic clouds, vapor cells or rare-earth doped 
crystals for quantum technology and -communication ap- 
plications. 
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